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Abstract : 

An analysis of MHD wave propagating in a gravitating and rotating medium perme- 
ated by non-uniform magnetic field has been done. It has been found that the Gradient of 
Magnetic Field when coupled with Rotation becomes capable to generate few instabilities 
(Temporal or Spatial) leading to the damping or amplification of MHD waves. The Jean's 
criterion is not sufficient for stability always. Rather, the waves will suffer instability unless 
their wave length (frequency) is less (greater) than certain critical values. Otherwise, those 
will smoothly propagate outward. Out of different scenarioes depending on the direction of 
the magnetic field, its gradient, rotation and wave propagation three important Special Cases 
have been discussed and different stability criteria have been derived. 

Finally, using the above theory we have obtained the stability/instability criteria for the 
waves moving parallel and perpendicular to the galactic plane in the Core and Periphery 
of the Central Region of Galaxy (C.R.G.) due to the coupled action of Rotation and 
Non-Uniform Magnetic field. The possibility of heating or occuring diffused condition 
inside the central region by MHD waves or smooth propagation of these waves (under some 
restrictions) through the C.R.G. has been briefly discussed. The numerical values of the pa- 
rameters of those waves for instabilities or smooth propagation have been estimated roughly. 
One may find some clues for the formation of Halo and Spiral Arms. 

keywords: MHD Instability, Magnetic field gradient and rotation, Central region of 
galaxy. 
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1. INTRODUCTION 



Any disturbance propagating through a non-uniform MHD fluid or plasma may be stable or 
undergoes instability. The damping of MHD waves is one of the causes of heating of such 
medium (Sturrock, 1966; Barnes, 1968; Bondyopadhaya et. al. 1972, 1974). In general, 
if streaming flow is present it may help instability in this medium (Sturrock, 1960; Gold, 
1965; Sudan, 1965). The strong rotation perpendicular to gravity, however, may help the 
stability of magnetized fluid (Gilman, 1970). The rotation can influence in the damping or 
heating of MHD medium like stellar interior (Bondyopadhaya 1972, 1974, 1978). The MHD 
instability due to temperature gradient and possibility of mass-outflow from the central region 
of galaxy has been investigated by Chakraborty and Bondyopadhaya (1998), Chakraborty 
(1998, Thesis). Sarkar and Bondyopadhaya (2007) have discussed the MHD instability due to 
non-uniform magnetic field and its effect on the propagation of waves of wave length greater 
than certain critical value. They also discussed the role of MHD instability in heating due 
to damping of waves having length less than that value. In reality, MHD flow seems to be 
connected with different physical processes of AGN (See, e.g. Blandford, 1990) or MHD Jets 
from AGN (See, e.g. Rosen and Hardee, 2000; Wiita, 2002). 

In this paper we have considered a gravitating and rotating MHD medium permeated 
by non-uniform magnetic field. We shall use the perturbation technique to linearise the 
guiding equations. Assuming the wave structure of the perturbation, General Dispersion 
Relation will be derived. Next investigation will be made about the wave propagation in 
different situations characterised by the Magnetic Field Gradient and Rotation. A number 
of stability/instability criteria will be derived which, specialy, are produced by the coupled 
action of Non-Uniform Magnetic Field and Rotation. Finally, these results will be used to 
have better understanding of instabilities (heating or diffused condition or other phenomena) 
or stabilities (smooth wave propagation) occuring in Galactic Central Region. 

2. BASIC EQUATIONS 

Let us consider a gravitating, rotating MHD fluid medium having finite conductivity. The 
basic equations of this fluid flow are the following : — - 



du/dt = 1/(4tt/xp)(V xB)xB-(Vxm)x u- p _1 Vp 

- V(u 2 /2 + (j)) -2(w x u) [Equation of Motion] (1) 

dB/dt = Vx(uxB) +r]V 2 B [MHD Field Equation] (2) 

dp/dt = —V(pu) [Continuity Equation] (3) 

\7 2 (f) = 47tGp [Poisson's Equation] (4) 

p = (C 2 /<y)p [Equation of State] (5) 



where, 

u = Fluid velocity 

B = Magnetic field 

p = Gas density 

<fi = Gravitational potential 

p = Permeability 

C s = Sound speed 

7=Ratio of specific heats 

All other symbols have there usual 



w = Angular velocity 

p = Hydrostatic pressure 

G = Gravitational constant 

rj = c/4irpa = Electrical Resistivity 

a = Electrical Conductivity 

c= Velocity of light in vaccum 



meaning. 
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3. PERTURBED LINEARISED EQUATION 

We shall now consider the situation when initially all the variables was at equilibrium but 
subsequently made perturbed as ip = ipo + ip/ where ipo and ip' are respectively unchanged 
value and small perturbation value of the variable ip . First all the variables of equations (1) 
to (5) are perturbed, then the equilibrium conditions (obtained by putting equilibrium value 
in (1) to (5)) are used and finally the square or higher power of perturbation variables are 
neglected so that we get the following linearised equations ( assuming, however, there is no 
initial streamflow i.e. uq = and 7 = 1) (See, Sarkar and Bondyopadhaya, 2007) 



du/dt = a[(V xB)x5 |(Vx B ) x B] - (p x )Vp - V</> - 2(w x u) (6) 

dB/dt = V x(uxB ) + 7]V 2 B (7) 

dp/dt = - Po V ■ u (8) 

V 2 </> = AirGp (9) 

p = clp (10) 



where 

C 2 = (RTq), To=Initial Temperature, a = l/(Airppo) 

R= Gas constant (obtained by dividing Universal Gas constant by mean mol. wt. of the gas). 
Here the primes have been droped and suffixes zero denote the equilibrium values. 

4. PERTURBATION STRUCTURE AND GUIDING EQUATION 

Let us assume that all the perturbation parameters are proportional to Exp i(kx-wt) i.e. they 
posses wave like structure , having the direction of propagation as x-axis, k and to being the 
wave number and wave frequency respectively. Then using the equations (8), (9) and (10) we 
can write the equations (6) and (7) componentwise as follows : 

From Equation of Motion : 

u x (uj 2 - k 2 c 2 + w 2 )/uj + u y (-2iw 0z ) + u z (2iw 0y ) + 

B y [a(-kB 0y +i(dB 0y /dx-dB Qx /dy))}+B z [a(-kB 0z -i(dB 0x /dz-dB 0z /dx))] = (11) 
u x (2iw 0z ) + u y {io) + u z (-2iw 0x ) + 

B x [-ia(dBoy/dx - dB 0x /dy)] + B y (kaB 0x ) + B z [ia(dB 0z /dy - dB 0y /dz)} = (12) 
u x (-2iw 0y ) + u y (2iw 0x ) + u z (u) + 

B x [ia(dB 0x /dz-dB 0z /dx)} + B y [-ia(dBo z /dy-dB 0y /dz)} + B z (kaB 0x ) = (13) 
From Field Equation : 

u x [-i(dB 0y /dy + dB 0z /dz)] + u y (idBo x /dy) + u z (idB 0x /dz) + B x (Lo + ir l k 2 ) = (14) 

u x (-kB 0y + idB 0y /dx) + u y [kB 0x - i(dB 0x /dx + dB 0z /dz)} + u z (idB 0y /dz) + 

B y (uj + i7]k 2 ) = (15) 

u x (-kB 0z + idB 0z /dx) + u y (idB 0z /dy) + u z [kB 0x - i(dB 0x /dx + dB 0y /dy)} + 

B z (L0 + ii]k 2 ) = (16) 

The equations can be put in matrix form as, 
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AX=0 (17) 

where, X=(u x u y u z B x B y B Z ) T is a column matrix and A = (a-ij) is a square ma- 
trix of order 6, denoting the coefficient of j-th variable in the i-th equation, i,j=l to 6, 
has the following values: 

a u = {lj 2 - k 2 c 2 s + w 2 )/u, a 12 = -2iw 0z , a 13 = 2iw 0y , a u = 0, 
ai5 = a[-kB 0y + i(dB 0y /dx - dB 0x /dy)], a 16 = a[-kB 0z - i(dB 0x /dz - dB 0z /dx)} 
a 2 \ = 2iw 0z , a 2 2 = u, a 2 3 = -2iw 0x , a 24 = -ia(dB 0y /dx - dB 0x /dy), a 25 = kaB 0x , 
a 2Q = ia(dB 0z /dy - dB 0y /dz) 

a 3 i = -2iw 0y , a 32 = 2iw 0x , a 33 = lo, o 34 = ia(dB 0x /dz - dB 0z /dx), 
035 = -ia(dB 0z /dy - dB 0y /dz),a 3 e = kaB 0x 

041 = -i(dB 0y /dy + dB 0z /dz), a 42 = idB 0x /dy, a 43 = idB 0x /dz, a 44 = u + irjk 2 , 
a 45 = 0, a 46 = 

a 5 i = -kB 0y + idB y/dx, a 52 = kB 0x - i(dB 0x /dx + dB 0z /dz), a 53 = idB Qy /dz, a 54 = 0, 
055 = a 44 , = 

a 6 i = -kB 0z + idB 0z /dx, a 62 = idB 0z /dy, a 63 = fc-B 0s - i(dB 0x /dx + dB Qy /dy), a 64 = 0, 
065 = 0, a66 = 055 = a 44 . 

where, = AirGpo = Gravitational potential. 

5. GENERAL DISPERSION RELATION (G.D.R.) 

Eliminating the variables from the equation (17) we get the General Dispersion Relation 
(G.D.R.) as, 

det(A)=0 (18) 

This is a relation of degree six representing six modes of wave propagation. 

6. CONDITIONS FOR INSTABILITIES 

We shall now analyse the wave instabilities under different conditions depending on direction 
of magnetic field. 

6.1 The magnetic field in the direction of wave propagation i.e. Bo = (Bo x , 0, 0). 
The dispersion relation (18) reduces to , 

J{A 2 -^ 2 w 2 x )-^A{w 2 y + w 2 z ) + JAV 2 x {L y 2 + L^ 
2uJu>kVj{ x [L y 1 (ujw 0z - 2iw 0x w 0y ) - L^^ujwoy + 2iw 0x w 0z )] - 

AuuV^woyL^+wozL- 1 ) 2 = (19) 

where k = k - i/L x , J = oj 2 - k 2 c 2 + wjj, V Ax = B 0x /(47r^po) 1/ ' 2 = Alfven velocity, 
A = lulu — kkV^ x ,L s being the characteristic length of variation of magnetic field 
along s i.e. x, y and z direction for s=x,y,z. It would be interesting to note that 2nd 
term shows that the waves will be affected by the joint action of non-uniformity of magnetic 
field along the wave propagation and rotation component perpendicular to wave propagation. 
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Furthar, the 4th, 5th and 6th term indicate that the waves which propagate along the mag- 
netic field but perpendicular to non-uniformity of magnetic field (i.e. L y , L z = finite) as well 
as rotation components (i.e wo y , w$ z / 0), must be affected by joint action of non-uniformity 
and rotation. Thus the joint action on the wave propagation along the magnetic 
field requires rotation component must be perpendicular to the wave propagation 
but non- uniformity may be parallel or perpendicular to wave propagation. 

Special case 1 : We assume that the waves are propagating along the gradient of mag- 
netic field (i.e. L x =finite but L y ,L z — ► oo ) and rotation component is in the direction of 
wave propagation only i.e. wq = (wq x ,0,0) (See Fig. 1) then from relation (19), 



J = 

uj{uj^2w Gx )/k 2 = Vl x [l-i/{kL x )\ 

Now relation (19.1) and (19.2) yield [taking n 

u? = k 2 C 2 - w 2 

k = ±[m — n 2 ] 1 / 2 + in 

where 

m = u(u =F 2wq x )/VI x , n = L~ l /2 



(19.1) 
(19.2) 

i.e to — ► u>] the following, 

(19.3) 
(19.4) 



From (19.3) one can observe that the temporal stable mode of waves propagating along 
the common direction of magnetic field, its gradient and rotation requires: 

Wave length , A < 2j ^f JL = A j( Jean's wave length) 

phase velocity, (%) = C 8 [1- ^j] 1/2 <C S = Sound Speed. (I a) 

The second relation (19.4) shows that the waves can avoid intsability provided its wave 
length A << 4ttL x 

Otherwise, such waves will undergo damping in the direction of increasing magnetic field 
[n > 0). This instability will be further enhanced at least for one mode whose 
wave frequency ( uj )< 2 rotational frequency (wq x ). In this situation the Instabilty Factor 
is given by- 

(3 = [± «(± 2 ™o*-") + (^) 2 ]V2 + l£ (lb) 

All these exhibit that the presence of magnetic field and rotation in the direction of 
wave propagation imposes restrictions on the stable and unstable mode's frequency and wave 
length. 

Special case 2 : In this case we assume the magnetic field acts along the direction of 
wave propagation, but its gradient and rotation act perpendicular to the wave propagation 
direction i.e. Bq = (Bo x , 0, 0), L x , L y — ► oo,wq = (0,0, wq z ) (See, Fig. 2) then relation (19) 
reduces to, 

uu + Vl x {L- 2 -k 2 ) = Q (20.1) 
J(wti - k 2 V\ x ) - Auju)wI z = (20.2) 
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We may note that the 1st relation involves non-uniformity of magnetic field but is independent 
of rotation while the 2nd relation involves rotation only but is independent of non-uniformity 
of magnetic field. Now, the relation (20.1) yields (taking i] — ► i.e. u — ► oS) , 

^ = V 2 x {k 2 - L- 2 ) 

Or, k 2 = (u> 2 + V 2 x L z 2 )/V 2 x >0 (20.3) 

These indicate that one mode will be temporally stable for 
A < 2ttL z . 

The other mode is governed by (20.2) which, in turn, can be written as, 

>2n2 , „.,2 a„.,2 7„2 T /2 1 , ,„2t/2 \u2r<2 „.,2l 



w 4 + ^\-k 2 C 2 s + w 2 - Aw 2 z - k 2 V% x ] + k 2 Vi x \k 2 C 2 - W 2 g ] = 



Or, k ^C 2 s Vl)-k 2 [u; 2 ^ + Vl(u J 2 + w 2 g )}+u; 2 (u J 2 + W 2 -4wl)=0 (20.4) 



If to 2 and ^2 are two roots of (20.4) then, 

uw 2 = evi x [k 2 c*-wi\ 



Evidently, two roots to 2 and to 2 are positive i.e the mode will be temporal stable if 
k 2 C 2 > w 2 i.e A < Aj= Jean's wave length (See, Ferraro and Plumpton, 1966). 

If k\ and k\ are two roots of (20.4) then, 
k\ + k 2 = [u 2 C 2 + V 2 Ax {u? + w 2 )]/(C 2 Vl) 
k 2 .k 2 = ^ 2 + w 2 g -A W 2 z )/{C 2 V 2 Ax ) 

Clearly the two roots k\ and k\ are positive i.e the mode will be spatially stable provided 

UJ 2 + W] > 4^4 

Thus the stable mode must have 

A < 2irL z , A < A j or uj 2 +w 2 > 4w^ z (II) 

Note that the second condition is Jean's condition. Thus the stable waves will be of higher 
frequency and of lesser wave length. Here both rotation and non-uniform magnetic field has 
been found to constrain the frequency and wave length but separately. 

6.2 Magnetic field perpendicular to wave propagation Here we will investigate 
the waves moving in the direction of the magnetic field but field gradient and rotation act 
perpenducular to the direction of wave propagation i.e. Bq = (0, Bq v , 0) or Bq = (0, 0, Bq z ) . 

6.2.1 For the first situation the dispersion relation (18) reduces to, 

J[uj 2 {lo 2 - 4w 2 x ) + uV 2 y L-\uL^ - ^I" 1 )] - 4ui 2 u 2 (w 2 y + <) 

+ 2ujkL y 1 L z 1 (Vl y ) 2 (w 0z L z 1 + ikw 0x ) - 4cvu}w 0z L z 1 V A : y (wo y L y 1 + w 0z L z 1 ) - 

ijju)V\ y {uo 2 -Awl x )k 2 '-\-2iuuL~ 1 L y l (ujw 0z +2iwo x w 0y )-8iuju;w 0x wozL~ 1 V Ay k = (21) 

where Va v = B$ y / (Ait pp®) 1 / 2 , L s being the characteristic length of variation of mag- 
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netic field along s i.e. x, y and z direction for s=x,y,z. 

It is interesting to note that there are product terms representing combined effect of rotation 
and magnetic field gradients (for example, 2nd part of 1st and 5th term, the 3rd, 4th, 6th 
and last term). 

If we analyse the relation (21) we can arrive at the conclusion that for gradient of magnetic 
field along the direction of wave propagation, the combined effect is peroduced 
provided wq x is present. However, the gradient of magnetic field along the direc- 
tion of magnetic field can not produce combined effect. Again if the gradient of 
magnetic field is in the z-direction, then the combined effect is present for wq z 
only. 

Special case 3 : Here we investigate the waves propagating (x-direction) perpendicular 
to magnetic field -Bo = (0,Bo y ,0) but its gradient and rotation both act perpendicular to 
both magnetic field and wave propagation (i.e. z-direction). Therefore, L x ,L y — > oo , 
w = (0,0, wqz) (See, Fig. 3) then the relation (21) reduces to, 

(J - 4wl)(uu + V 2 y L~ 2 ) = io 2 k 2 V\ y (21.1) 

This relation reveals that the combined effect could be produced by magnetic field 
gradient and rotation provided both act perpendicular to magnetic field as well 
as wave propagaton. 

Now, relation (21.1) yields (taking r\ — ► 0), 

w 4 + w 2 [Vl y {L~ 2 - k 2 ) - (k 2 C 2 s -w 2 g + Awl)} + Vl y L-\-k 2 Cl + w 2 - Aw 2 z ] = 
Or k 2 = i^+V 2 L 7 2 )(u 2 +w 2 g -Awl)/[co 2 V 2 +C 2 (u J 2 +V 2 L z 2 )} (21.2) 



If uj\ and u>2 are two roots of the first relation in (21.2) then, 
<"? + <4 = Vi y (k 2 - L- 2 ) + k 2 C 2 - w 2 + Awl 

= k 2 (C 2 + V 2 y ) - {V 2 Ay L- 2 + w 2 g - Awl) 
uW 2 = -V 2 L-\k 2 C 2 -wl + Awl) 



It is to be noted that if k 2 C 2 > w 2 i.e. A < Xj, then both the roots can not be posi- 
tive. This Jean's criteria can not ensure stable mode so long as magnetic field, its gradient 
and wave propagation are mutually perpendicular to each other (See, Fig. 3). However, if the 
non-uniformity of magnetic field disappears i.e. L z — > oo, then the first relation of (21.2) 
yields, 

u J 2 = k 2 (V 2 y + C 2 )-w 2 g +Aw 2 z 

This shows that for temporal stability Jean's criteron is sufficient but for spatial stability 
the requirement is to 2 + w 2 > Awl 

Clearly, the temporal stability of the waves require that both the roots (uj 2 and oj^) must be 
positive. This means 

k 2 (C 2 + V 2 y )>(V 2 y L z 2 + w 2 -Aw 2 z ) i.e A < 2n[ y2 °l + J}_ Au]2 Y' 2 

Ay 2 9 Oz 

and 
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k 2 C 2 -w 2 + Aw 2 z < i.e A > 2n [wj _^2 J1/2 

These give a range of wave length for the temporal stability of waves. Namely, 

2tt[ - 2 T °t +V f - 2 ] l l 2 > A > 2tt , 2 fta 11/2 (III) 

Evidently, the rotation as well as magnetic field gradient both are capable to impose re- 
strictions on the upper bound of the wave length of temporally stable modes while only the 
rotation is capable to restrict the lower bound of the temporally stable mode. 

6.2.2 For the second situation, the dispersion relation (18) reduces to , 

J[uj\uj 2 - 4w 2 x ) + uV 2 z L-i(uL-i + 2iw 0x L~ 1 )} - ^\w 2 y + <) 

- 2u}kLy l L- l (Vl z ) 2 (w 0y Ly 1 + ikw 0x ) - ^LOWoyLy X V\ z {woyLy X + wozLj 1 ) - 

loloV^ z (lo 2 —4wQ X )k 2 '^iLULuL' 1 L^ 1 (iJWoy—2iwo x wo z )—8iuJujwo x woyLy l Vj z k = (21-3) 

where Va z = Bqz/^tthpo) 1 / 2 , L s being the characteristic length of variation of mag- 
netic field along s i.e. x, y and z direction for s=x,y,z. 

The conclusions which can be drawn from the relation (21.3)are similar to that from (21). 

SUMMERY : We have discussed a few instabilities of the many cases which could be 
covered by the dispersion relation (18). For such cases pertaining to different situations 
we can study the instabilities, the role of combined effect of rotation and non-uniformity of 
magnetic field. But these will be done elsewhere. At present we can summerise the above 
study as follows. 

When the wave propagates along x-direction we have considered three cases results of 
which are the following. 



Special 

Cases 


B 


VB 


w 


Dispersion Relation 
* Stability Condition 


Figures 


1 


Bq x 




WOx 


J = 

w(w T 2w 0x )/k 2 = V\ x [\ - i/(kL x )} 
* A < Aj or A << 4irL x 


Fig. 1 


2 


Box 






louj + VI x {L~ 2 -k 2 ) = 
J{ujuj - k 2 V\ x ) - Auouowl z = 

* oj 2 + w 2 > 4wq z 
* or A < 2-kLz or A < Aj 


Fig. 2 


3 


B y 


V ' z Bo y 


WQz 


(j - 40(^ + VlyLT 2 ) = ^k^y 

* oj 2 + w 2 > 4wq z 
* or 2tt[— 2 T °X +V f - a Y' 2 > A > 2tt r 2 f° 2 11/2 


Fig. 3 



It is found that Jean's criterion is necessary but not sufficient for the temporal stability of 
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waves in Special Cases 1 and 2. However, the waves with shorter length and higher frequency 
are generally stable. For the Spatial stability the condition for Cases 2 and 3 are same. The 
Temporal stability for case 3 is too restricted to exist so long as non-uniformity of magnetic 
field exits. 

7. CENTRAL REGION OF GALAXIES 

The theory which has been discussed so far may be useful to explain different astrophysical 
phenomena like heating of the medium and mass outflow from it, provided, the medium could 
be treated essentially as MHD or Plasma. 

Let us first note down the characteristic features of the medium embeded in the central region 
of galaxies (See, e.g. Chakraborty and Bondyopadhaya,1998). 

1. Fluid characteristics : 

We know that the ionised gas medium could be treated as fluid if the mean free path of the 
charged particles is less than the linear dimension of the medium ie. A m < L (See, e.g. Wolt- 
jer, 1965). Now observations have revealed that huge material in different galactic central 
regions exists which may satisfy this criterion (See, e.g. Krotik et al., 1981; Perola et al., 
1986; Uberoi, 1988). We shall study the dynamics of this fluid in the central regions by the 
help of hydromagnetic theory developed in the text (See also, Chakraborty and Bondyopad- 
haya, 1998). 

2. Magnetic field: 

Woltjer (1965, 1971, 1990) discussed the important role of magnetic field in the activities 
of Galactic Nuclei, Spiral Arms, Seyfert Galaxies and Radio jets (See also, Zeilik and Niley, 
1997). This view was supported by many others (e.g. Osterbrock and Mathews, 1986). In- 
fact, the existence of ~ 1 gauss magnetic field in AGN has been reported by many authors 
(See, e.g. Rees, 1987, Wielebinski, 1990, Morris, 1990). There is, however, indication 
that the field is generally parallel to the plane in denser region (core) and per- 
pendicular to the plane in low density region (periphery) (See, Chuss et al., 2003). 
The magnetic field parallel or perpendicular to the galactic plane is likely to be non-uniform. 

3. Rotation: 

Generally rotation axis is either parallel or slightly inclined to the normal of the galactic 
plane. Therefore, for convenience we take the dominant component of rotation as perpendic- 
ular to the plane. 

4. Division of Central Region: 

There are two zones in the central region of the galaxy. 

i) Core Region: Obviously the core of the central region is a denser medium. Following 
Chuss et al., 2003, let us assume that this region has magnetic field parallel to galactic plane. 
As per discussion above the rotation can be taken as perpendicular to this galactic plane. 
Furthur, the magnetic field gradient can be taken parallel and perpendicular to the galactic 
plane. 

ii) Peripheral Region: Similarly, the peripheral region of the centre is a lighter medium. 
As per suggession Chuss et al., 2003, the magnetic field can be taken as perpendicular to the 
galactic plane. The magnetic field gradient and rotation can also be taken as before. 
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8. STABILITY OR INSTABILITY OF THE CENTRAL REGION AND 

SIMILAR MEDIA 



Now, we shall study the stability or instabilitiy of the wave propagation in the Core and 
Periphery of central region of the Galaxy. 

i) The Core region: It is of higher density where the magnetic field may be taken as 
parallel but its gradient may be taken as perpendicular to the galactic plane. The MHD 
wave propagation in such system has been discussed in Special Case 3. The combined effect 
of rotation and magnetic field gradient has been found on the waves propagating parallel to 
the galactic plane. In this case, it is observed that the waves will be spatially stable for 

4w 2 z <u 2 + w 2 . 

And, the wave length of temporally stable wave requires the constrain (III) to be satisfied. 
We may note that if L.H.S.is greater than R.H.S. then 

w 2 g > 4w 2 0z + C 2 /L 2 Z i.e. the density p > &3+^oX!l 
The condition (III) is true when the galactic plane flow (x-direction) is perpendicular to 
magnetic field (y-direction) . However, the flow along the magnetic field could be described 
by the theory outlined in Special Case 2 (Fig. 2). For stability such wave requires 

uj 2 + w 2 > 4wq z (as before) and A < Aj, A < 2nL z 

When all these stable waves parallel to the galactic plane, and both along and perpen- 
dicular to the magnetic field continue to propagate for million of years those may feed the 
Spiral Arms or even may lead to the formation of Spiral Arms or other outer structures of 
the galaxies including our Galaxy. 

ii) Peripheral Region : This region is of low density where the magnetic field, its 
gradient and rotation all act perpendicular to the galactic plane (as shown in Special Case 1, 
Fig. 1). Here, the waves propagating perpendicular to the galactic plane undergo instability 
due to the joint action of rotation and magnetic field gradient. One temporal stable mode 
can propagate whose phase velocity is less than sound velocity and whose wavelength A < Aj. 
The other mode may be stable provided A << 4ttL x . If the wavelength does not satisfy this 
condition the waves will be amplified in the direction of decreasing magnetic field . In fact, the 
magnetic field decreases away from the galactic plane. In addition the waves with frequency 
co < 2wq x will undergo instability with instability factor (I b) (as shown in Special Case 
1), otherwise, I.F.=1/(2L X ) for the waves with period r > tt/wq x these unstable waves can 
create a difnused condition of material similar to Halo. While the stable waves can propagate 
in upward direction. However, it is interesting to note that if there is any flow towards the 
galactic plane in the direction of increasing magnetic field, the waves will be damped leading 
to the heating of the medium. In the long run these may help to create hot Halo like Solar 
Corona. 

Now it is reported that from the galactic central region most of the energy come out in 
the form of infrared emission. In general, active galaxies have not only radio, u-v and x-ray 
but also infrared emission. In fact, later emission may be from heated dust or hot ionised 
gases (See, e.g. Zeilik, 1997). Heated dust around the nucleus of galaxies may absorb high 
energy radiation emitted from energised Core and re-emits it at larger infrared wave length 
(See, e.g. Moche, 2004). 
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Infact, any instability may initiate heating of the medium via damping of MHD waves. 
Of course, not all waves are allowed to do it. Generally, waves with shorter wave lengths can 
pass through and longer waves are damped or amplified ( See e.g. Sarkar and Bondyopad- 
haya, 2007). Again, the system rotating with very high velocity may cause the instability 
of the MHD waves (See, Bondyopadhaya, 1972). Therefore, both rotation and non-uniform 
magnetic field may be held responsible for the heating of the Central Region. But both in the 
Core and Periphery of the central region the waves with shorter wavelength or high frequency 
can move smoothly. For such waves moving parallel or perpendicular to the galactic plane, 
when continues for long time may be responsible for the formation of the Spiral Arms or Halo 
region respectively. 

We have seen in the text that the non-uniform magnetic field when coupled with rotation 
is capable to produce instability. This type of instability could not be produced by them 
separately (See, e.g. Barnes, 1967; Bondyopadhaya, 1972; Paul and Bondyopadhaya, 1973; 
Yusef-Zadeh et al., 1986; Chakraborty and Bondyopadhaya, 1998). That is to say, this 
instability is a product function of non-uniform magnetic field and rotation. In reality, many 
astrophysical media posses both rotation and non-uniform magnetic field (the rotation axis, 
however, may not coincide with magnetic axis always). Therefore, the theory like above may 
have good relevance for such media like AGN, Galactic Central Regions, Seyfert Galaxies 
and Infrared Galaxies. Hence, the joint effort of rotation and non-uniform magnetic field for 
initiating instability can be found there. In other words, the events occuring there may be 
at least partially due to joint action of non-uniform magnetic field and rotation. 

9. NUMERICAL ESTIMATION 

In the previous section we have discussed the instability criteria in the central region of 
galaxy (C.R.G.). Now, we will make a few numerical estimations of those criteria. We have 
considered the central region of our galaxy as a hot plasma bed described by the MHD equa- 
tions (1 to 5) of Sec. 2 (See e.g., Chakraborty and Bondyopadhaya, 1993, 1998; Sarkar and 
Bondyopadhaya, 2007). Let us first consider some physical parameters in the central region 
of our Galaxy. 

The radius of the central region may be taken as 0.5 to 100 pc because well defined cir- 
cumnuclear gas disks have been found upto that region in a number of galaxies (Morris, 1998). 

Further, the central region may be considered as a medium of uniform density, which is 
given by po ~ 3.5 x 10~ 20 g cm -3 i.e. 10 36 g pc~ 3 (See Balick and Sanders, 1974; Genzel et 
al., 1985; Lo, 1986 ). Consequently, gravitational frequency becomes w g = (AirGpo) 1 / 2 ps 
1.7 x 10~ 13 s -1 ps 5.36 per mil. yr. 

The adiabatic sound speed is given by C s ss Q.5pc per mil yr (for To = 10 4 i^ at 10 pc 
from the centre, vide Yusef-Zadeh et al., 1984). 

The strength of general magnetic field is 4-6 p,G with a cell size of 10-100 pc (See, Ohno 
and Shibata, 1993). Let us take magnetic field strength as 5pG and the variation of magnetic 
field as 2/j.G over 50 pc range. Thus the characteristic length of variation L x or L z ss 25pc. 
Hence, Alfven velocity Va ps 1.4 x 10~ 18 pc s -1 ss 4.41 x 10~ 5 pc per mil. yr . 
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The central region of the Galaxy rotates like a rigid body. Beyond this region the ve- 
locity first decreases and then increases. At nearly 8 kpc from the centre the velocity be- 
comes maximum. The Sun is situated at 8.5 kpc from the centre having linear velocity 
220fcms~ 1 . Then Sun's angular velocity about galactic centre ~ 84.6 x 10 17 s . (See, Kar- 
tunnen, 2007). Let us take this value as same as the angular velocity of central region i.e. 
w 0z ps 84.6 x lO-^s- 1 xi 2.7 x 10" 2 per mil. yr. 

Therefore, we obtain a critical wave length (i.e. Jeans wave length) 
Aj = 2TTC s /wg 7.62pc. 

Next, let us find what are the waves which can propagate or which will suffer instability 
inside the two areas of the central region of our Galaxy. 

i) Peripheral Region : (Special Case 1, Fig. 1) 

Here the magnetic field, its gradient and rotation all act perpendicular to the galactic 
plane. One modes of the waves moving perpendicular to this plane will be temporal stable 
provided 

A < Aj ~7.6 pc and the other mode require 

A << 4irL x ~ 315 pc. Conversely if the waves have length greater than 7.6 pc those 
will suffer instability. In addition, if A is not much less than 315 pc some more modes will 
be unstable. Moreover, the waves having frequency less than 2wq x ~ 17 x 10 _16 s _1 (i.e. 
r > 10 8 yr) the spatial instability will be enhanced. The actual instability factor is given by 
f3 in Special case 1. 

Such waves will be amplified leading to the defused condition if it moves 
along the decreasing direction of magnetic field i.e. away from the Galactic plane 
towords Halo. If, however, there is inflow of MHD waves from outside towards 
the Galactic plane i.e. in the direction of increasing magnetic field then the waves 
will be damped leading to the heating of the Halo . 

ii) Core Region : (Special case 2,3. Fig 2,3) 

The spatial stability of waves in this region require w g > 2wq z . As we have already 
discussed here w g ~ 1.7 x 10~ 13 s -1 and wq z ~ 8.5 x 10 6 s . Thus the inequality is 
satisfied. Hence, waves of all frequencies will spatially stable. 

Now, the waves will propagate, in the direction perpendicular to the magnetic field, 
smoothly (i.e. temporally stable) must have the wave length restricted by the following 
inequality 

2tt\ C ° +V Ay 11/2 > x > 2tt ^ 

Incidentally, w g (x 1.7 x 10~ 13 ) >> 2w 0z (x. 1.7 x 10~ 15 ) and 
C s » Va v , Va v /L z «3x 10 _20 s _1 . These values will make R.H.S and L.H.S of the above 
inequality almost same. Therefore, no temporally stable wave can exist. However, those 
waves moving along the magnetic field will be temporally stable whose wave length is given 
by 



3 We make a note of the fact that the source of energy of the solar corona is still not known but it is thought 
to be some how associated with magnetic reconnection and, therefore must be understood in terms of MHD 
process (See, Cravens, 1997). 
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A < 2ttL z 157. 5pc. 



As we have seen that in the central region of our Galaxy w g » 2wq z and 
C s » Va v but the situation might be different in the long past in the life of our 
Galaxy itself. When both magnetic field and rotation was much higher. Similar 
situation may be present in some other galaxies also. For example, in AGN 
magnetic field is 10 6 times (See, Sec. 7) the field that we have taken, so that the 
contribution due to magnetic field or its gradient is much more significant. The 
derived condition could be applicable in those galaxies particularly at the time 
of formation. 

10. CONCLUSION 

In this paper the General Dispersion Relation (G.D.R.) for the unidirectional (x-direction) 
wave propagation has been derived in the rotating medium where magnetic field is non- 
uniform. Then the guiding dispersion relations have also been obtained for magnetic field 
either along or perpendicular to wave propagation. There could be many different situations 
depending on the magnetic field, its gradient and the direction of rotation. But out of these 
only three Special Cases have been discussed. 

A number of temporal/spatial instability /stability conditions have been derived. Com- 
paring the Special Cases it is observed that the combined effect of rotation and mag- 
netic field gradient is present in the product form provided the wave propagates 
along the magnetic field. The coupled action of rotation and magntic field's non- 
uniformity have also been noticed when the wave propagates perpendicular to 
the magnetic field and its gradient. In general the instability conditions reveal that the 
waves having shorter frequency or longer wave length become unstable. 

The instability conditions, thus derived, could be used to study different phenomena 
occuring in the the medium similar to the central region of our Galaxy. It is suggested that 
such stabilities or instabilities are more significant in the early stage of evolution of many 
galaxies when both rotation and magnetic field are much higher. Those waves which can 
propagate from the central region may help in the formation of Spiral Arms and Halo of the 
galaxy in long run. 

REMARKS 

i) Mathematical analysis have been given more stress in this paper. 

ii) A number of sub cases have also been analysed but all those are not included in this 
paper. 

iii) The theory of the formation of hot Halo and Spiral Arms from the Galactic central 
material appears to be consistence with the underlined proposition of this paper. 

iv) It is well known that the electrical resistivity causes the wave unstable. If r\ does not tend 
to 0, the dispersion relations obtained in the text could be easily derived. This will, however, 
be considered in separate communication. 
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